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Abstract 



Using simple space-time implementation of the random cascade model we investigate nu- 
merically influence of the possible fractal structure of the emitting source on Bose-Einstein 
l/~) . correlations between identical particles. The results are then discussed in terms of the non- 

extensive Tsallis statistics. 
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1.1 Introduction 

Two features seen in the analysis of multiparticle spectra of secondaries produced in high energy 
collision processes are of particular interest: (?) intermittent behaviour observed in analysis of 
factorial moments and (u) Bose-Einstein correlations (BEC) observed between identical particles. 
Whereas the former indicates a possible (multi) fractal structure of the production process (in 
the momentum space) jj] the latter provides us with knowledge on the space-time aspects of 
production processes 0] . It was argued [H that these features are compatible with each other only 
when: (i) either the shape of interaction region is regular but its size fluctuates from event to event 
according to some power-like scaling law or (ii) the interaction region itself is a self-similar fractal 
(in coordinate space) extending over a very large volume. Although there exists a vast literature 
on the (multi) fractality in momentum space B its space-time aspects are not yet fully recognized 
with |3J remaining so far the only representative investigations in this field. We shall present here 
numerical analysis of a particle production model possessing both the momentum and coordinate 
space fractalities g extending therefore ideas discussed in || to a more realistic scenario. 
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1.2 Cascade model used 

As a model we shall choose a simple self-similar cascade process of the type discussed in Q but de- 
veloped further to meet our demands. The following point must be stressed. Every cascade model 
is expected to lead automatically to intermittent behaviour of momentum spectra of observed par- 
ticles J7J. Although this is true for models based on random multiplicative processes in observed 
variables (like energy, rapidity or azimuthal angle), this is not necessarily the case for multiplica- 
tive processes in variables which are not directly mesurable but which are, nevertheless, of great 
dynamical importance (like masses Mi of intermediate objects in a cascade process considered 
here). In a purely mathematical case, where cascade process proceeds ad infinitum, one eventually 
always arrives at some fractal picture of the production process. However, both the finite masses fj, 
of produced secondaries and the limited energy M stored originally in the emitting source prevent 
the full development of such fractal structure j8| . One must be therefore satisfied with only limited 
and indirect presence of such structure. This applies also to the analysis presented here. 

In our model some initial mass M "decays" into two masses, M — > Mi + M 2 with M 12 = 
fci.2 • M and ki + k 2 < 1 (i.e., a part of M equal to (1 — k\ — k 2 )M is transformed to kinetic energies 
of the decay products M\ t2 ). The process repeats itself until Mi. 2 > M [}X being the mass of the 
produced particles) with successive branchings occuring sequentially and independently from each 
other, and with a priori different values of ki,2 at each branching, but with energy-momentum 
conservation imposed at each step. For different choices of dimensionality D of cascade process, 
D = 1 (linear) or D = 3 (isotropic), and for different (mostly random) choices of decay parameters 
k\ t 2 at each vertex, we are covering a variety of different possible production schemes, ranging from 
one-dimensional strings to isotropic three-dimensional thermal-like fireballs. For our purpose of 
investigation of connections between BEC and space-time fractality of the source we have extended 
this (momentum space) cascade also to the space-time and we have added to it a kind of BEC 
"afterburner" along the lines advocated recently in M. 

What concerns space-time development we model it by introducing a fictitious finite life time 
t for each vertex mass Mi, distributed according to some prescribed distribution law given by 
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This procedure is purely classical, i.e., intermediate masses Mi are not treated as resonances (as 
was done in CO]) but are regarded to be stable clusters with masses given by the corresponding 
values of decay parameters k\^ 2 and with velocities (3 — P\, 2 jE\^ 2 ({Ei,2', Pi, 2) are the corresponding 
energy-momenta of decay products calculated in each vertex in the rest frame of the parent mass) . 
The energy-momentum and charges are strictly conserved in each vertex separately. The form of 
T used in (|l|) allows to account for the possible fluctuations of the evolution parameter r [[llj and 
finds its justification in the Tsallis statistics, to which we shall return later when discussing our 
results @. 

1.3 Bose-Einstein correlations 

Our main goal is the investigation of the BEC, in particular whether these correlations show 
indeed some special features which could be attributed to the branchings and to their space-time 
and momentum space structure. We are therefore interested in two particle correlation function 



To calculate it we have decided to use the ideas of the BEC "afterburners" advocated recently in 
U . Such step is necessary because cascade per se do not show bosonic bunching in momenta (as is 
the case in models where Bose statistics is incorporated from the beginning, like |Q3; however, we 
cannot follow this strategy here). Because we are interested in some possible systematics of results 
rather then in particular values of the "radius" R and "coherence" A parameters characterizing 
source M, we have chosen the simplest, classical version of such afterburner. After generating a 
set of i = 1, . . . , Ni particles for the I th event we choose all pairs of the same sign and endow them 
with the weight factors of the form 

C = 1 + cos [(n - r 3 ){jpi-pj)] (3) 

where r, = (ti,fi) and pi — (Ei,pi) for a given particle. The signs are connected with charges 
which each cascade vertex is endowed with using simple rules: {0} — ► {+} + { — }, {+} — > {+} + {0} 
and{-}^{0} + {-}. 

2 Results 

Although it is straightforward to get our cascade model in the form of the Monte Carlo code, the 
main features of D — 1 case can be also demonstrated analytically. For example, in the limiting 
cases of totally symmetric cascades (where for all vertices k\ t 2 = k), in which amount of energy 
allocated to the production is maximal, one gets the following multiplicity of produced particles: 



It is entirely given by the length of the cascade, L max = ln(M / /i) / \n(l / k) , /i = \rrr% + (pt) 2 - 
The exponent d,F is formaly nothing but a generalized (fractal) dimension of the fractal structure 
in phase space formed by our cascade. Notice the characteristic power-like behaviour of N S (M) in 
(0) which is normally atributed to thermal models. For example, for k = 1/4 one has N s ~ A/ 1 / 2 , 
which in thermal models would correspond to the ideal gas equation of state with velocity of sound 
Co = l/v3 |5[|. In the opposite limiting case of maximally asymmetric cascades, M — > /i + Mi 
(where fci = fx/M and ki — k) in which the amount of kinetic energy allocated to the produced 
secondaries is maximal, the corresponding multiplicity is equal 

1 M 

N a = 1 + L max = 1 + -— . In— . (5) 

ln i ^ 

The dependence on L max is now linear (i.e., dependence on the energy is logarithmic). The 
important feature, which turns out to be valid also in general, is the observed scaling in the ratio 
of the available mass of the source M and the mass of produced secondaries: M//J,. 

The D = 3 case differs only in that the decay products can flow in all possible directions, which 
are chosen randomly from the isotropic angular distribution. To allow for some nonzero transverse 
momentum in D = 1, we are using the transverse mass fi = 0.3 GeV. For D = 3 cascade fj, is 
instead put equal to the pion mass, \i = 0.14 GeV. All decays are described in the rest frame of the 
corresponding parent mass Mi in a given vertex. To get the final distributions one has to perform 
a necessary number of Lorentz transformations to the rest frame of the initial source mass M. As 
an output we are getting in each run (event) a number Nj of secondaries of mass fi with energy- 
momenta {Ej]Pj)i and birth space-time coordinates {tjjfj)., i = l,...,Nj (i.e., coordinates of 



the last branching). Results presented in Figs. 1 and 2 are obtained from 50000 such events. 
Decay parameters k\^ were chosen randomly from a triangle distribution P(k) = (1 — k) (leading 
to a commonly accepted energy behaviour of N(M) ~ ^0.4-^0.5^ as discussed above). For more 
detailed presentation of rapidity and multiplicity distributions and demonstration of intermittent 
behaviour of factorial moments see f| . 

Fig. 1 shows densities p(r) of points of production and correlation function C 2 (Q) as defined 
in Q3J) for D = 1 and 3 dimensional cascades originated from masses M = 10, 40 and 100 GeV. 
The evolution parameter is set equal r = 0.2 fm (in B we discuss also r ~ \/M case). The decay 
function T is taken exponential (i.e., q = 1). We observe, as expected in [0, a power-like behaviour 
of cascading source: 
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p(r) ~ I - J r > r , (6) 

but only for r > tq, i.e., for radii larger then some (not sharply defined) radius ro, value of which 
depends on all parameters used: mass M of the source, dimensionality D and evolution parameter 
t of the cascade. Below tq the p(r) is considerably bended, remaining almost flat for D = 1 
cascades. For the limiting case of M = 100 GeV the corresponding values of parameter L vary 
from L = 1.89 for one D — 1 cascades to L = 2.78 for D = 3 cascades. The shapes of p(r) scale in 
the ratio M/p in the same way as the multiplicities discussed before (the same remains also true 
for rapidity and multiplicity distributions, cf. ||). One can summarize this part by saying that 
power-like behaviour || sets in (albeit only approximately) only for long cascades (large values 
of M/p). It remains therefore to be checked whether (and to what extend) such conditions are 
indeed met in the usual hadronic processes. 

The corresponding BEC functions C2 show a substantial differences between D = 1 and D = 3 
dimensional cascades, both in their widths and shapes. Whereas the former are more exponential- 
like (except, perhaps, for small masses M) the latter are more gaussian-like with a noticeably 
tendence to flattening out at very small values of Q observed for small masses M. Also values 
of intercepts, C 2 {Q — 0), are noticeable lower for D = 3 cascades. The length of the cascade 
(i.e., the radius of the production region, cf. discussion of density p before) dictates the width 
of C%{Q). However, the M/p scaling observed before in shapes of source functions is lost here. 
This is because Ci depends on the differences of the momenta p = p coshy, which do not scale 
in M/p. The flattening mentioned above for D = 3 cascades are the most distinctive signature 
of the fractal structure combined with D = 3 dimensionality of the cascade. The correlations of 
the position-momentum type existing here as in all flow phenomena are, in the case of D = 3 
cascades, not necessarily vanishing for very small differences in positions or momenta between 
particles under consideration. The reason is that our space-time structure of the process can have 
in D = 3 a kind of "holes", i.e., regions in which the number of produced particles is very small. 
This is perhaps the most characteristic observation for fractal (i.e., cascade) processes of the type 
considered here. This feature seems to be more pronounced for diluted cascades corresponding to 
q < 1 case discussed below. 

Fig. 2 displays the same quantities but this time calculated for M = 40 GeV and for different 
values of parameter q in the decay function T defined in (|l|) . This function is written there in form 
of Tsallis distribution |l^] , which allows to account for a variety of possible influences caused by, 
for example, long-range correlations, memory effects or for the possible additional fractal structure 
present in the production process. They all result in a non-extensitivity of some thermodynamical 
variables (like entropy) with |1 — <?| being the measure of this non-extensitivity. In practical terms 



of interest here, for q < 1 the tail of distribution (|T|) is depleted and its range is limited to 
t S (0, t/(1 — q)) whereas for q > 1 it is enhanced in respect to the standard exponential decay law 
(and its range is t € (0, oo)) |L2j. In other words, one can account in this way for both more diluted 
(for q < 1) and more condensed (for q > 1) space-time structure of the developing cascades. Such 
distributions are ubiquitous in numerous phenomena and they are founded in the, so called, Tsallis 
non-extensive thermostatistics |12J generalizing the conventional Boltzmann-Gibbs one (which in 
this notation corresponds to the q = 1 case). It has found also applications in high energy and 
nuclear physics (cf. pi for references). Its effect on the cascades investigated here is, as can be 
seen from Fig. 2, in that it mimics (to some extent) the changes atributed in Fig. 1 to different 
energies (making cascade effectively shorter for q = 0.8 and longer for q = 1.2). The results of Fig. 

2 (taken for M — 40 GeV) should be then compared with those of Fig. 1 for M — 10 and M = 100 
GeV. They demonstrate that effects of longer or shorter cascades in momentum space (as given by 
different M in Fig. 1) is similar to effects of the more or less condensed cascades in the position 
space as given by q here. This fact should be always kept in mind in such analysis as ours. 

3 Conclusions 

We conclude that BEC are, indeed, substantially influenced by the fact that the production process 
is of the cascade type (both in momentum and space-time) as was anticipated in 0], although 
probably not to the extent expected (which, however, has not been quantified there). In practical 
applications (fitting of experimental data) there are many points which need further clarification. 
The most important is the fact that data are usually collected for a range of masses M and among 
directly produced particles are also resonances. This will directly affect lengths of the cascades, 
and through them the final results for Ci. Selecting events with similar masses M should allow 
to check whether in such processes q = 1 or not. The importance of this finding is in that q < 1 
would signal a non-stochastic development of the cascade, whereas q > 1 would indicate that, as 
has been discussed in [ |ll| , parameter r fluctuates with relative variance 

w = —i =q-l. 7 

Such fluctuations are changing exponential behaviour of V in (Ml) to a power-like distribution with 
enhanced tail. 
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Figure Captions 

Fig. 1. Density distribution of the production points p(r) (left panels) and the corresponding 
C2(Q = \Pi ~ Pj\) (right panels) for one-dimensional (upper panels) and three-dimensional (lower 
panels) cascades. Each panel shows results for M = 10, 40 and 100 GeV masses of the source. 
Time evolution parameter is r = 0.2 fm and nonextensitivity parameter q = 1. 
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Fig. 2. The same as in Fig. 1 except that each panel shows results for mass of the source M = 40 
GeV and for three different values of the nonextensitivity parameter q = 0.8, 1.0 and 1.2. 
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